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Consideration of the adiabatic character of radiation expansion in early Universe 
leads to the conclusion that equilibrium distribution of the primordial radiation in 
the presence of charged particles could be different from the Planck distribution 
in some regions of the spectrum. The equilibrium distribution of electromagnetic 
radiation (the black body radiation) is generalized for the system containing an 
extremely dense fully ionized plasma. The conditions of the adiabatic expansion of 
radiation for the model of the early Universe are found. 
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The equilibrium distribution of photons in the background cosmic radiation [1, 2] obeys 
the Planck law with high accuracy in the microwave region of the spectrum. However, the 
primordial radiation existed in the high density plasma medium. The plasma composition in 
early Universe changed with time. Evolution of the composition of the elementary particles 
in this plasma is not yet clear enough. However, it is possible to consider the changes of the 
temperature and density in the adiabatically expanding plasma of a fixed composition. 

The presence of the fully ionized plasma leads to the distortion of the Planck distribution 
of radiation in low- frequency limit, when the plasma frequency Q p is of the same order with 
the current frequency. The thermodynamic functions of such radiation are the functions of 
the dimensionless parameter hfl p /T. When hQ p /T <C 1, the thermodynamic functions for 
radiation closely coincide with the ones obtained with the Planck distribution. However, 
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in the case when hfl p /T > 1 the situation dramatically changes. Realization of this latter 
condition in the early Universe seems possible. In this case our understanding of the history 
of the Universe is modified. 

According to the Planck law [3, 4] the spectral energy distribution of the black radiation 
in volume V is equal 

V hu 3 du 

dE{uj) = ^exp(huj/T)-V (1) 

where c is the velocity of light. This formula implies that there is equilibrium between 
photons and matter at temperature T. Interaction between the photons and matter should 
be small enough as not to disturb ideality of the photon gas. The net absorption and 
creation of photons by matter in the volume V are absent, these processes are in a balance. 
Also, interaction between the photons themselves is extremely weak. At the same time, to 
reach equilibrium in the photon gas, a weak photon-matter interaction is necessary. These 
conditions are satisfied quite well in gases for all frequencies that lie far away from the 
frequencies of the absorption lines of matter. 

The Planck law has been verified experimentally, e.g. by pumping laser radiation into 
the empty (or filled with a rarefied gas) cavity through a small hole and then observing the 
spectral distribution of the radiation emitted through the hole. The spectral distribution of 
the emitted radiation is the Planck distribution with high accuracy. 

At the same time there exists a medium, plasma, in which the electromagnetic radiation 
can exist without damping or with a weak damping, and it possesses a different dispersion 
than that in vacuum, where uj = ck. 

Let us consider a fully ionized plasma in some volume V. We suppose that the plasma 
itself is in the equilibrium state, and having the same temperature T as the surrounding 
matter. The condition of small damping of the transversal mode uj = y / c J P~+~n| in plasma 
[5] is uj ^> v, where Vt p is the plasma frequency of the electrons and v is the characteristic 
collision frequency between the charged particles. The characteristic values of the collision 
frequency are equal to v — n 3 / 2 Ze^nL / A\f2T 3 l 2 m l l 2 for low frequencies, and to v u = 
4V27rZe 4 nL/3T 3 / 2 m 1//2 - for high frequencies. Here, v ~ where n is the electron density 
and L is the Coulomb logarithm for the electron-ion collisions. Influence of the collisions 
on the radiation field is small if Q p ^> and uj ^> v w . The condition of weak interaction 
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between the particles can then be written as 

3 

r « 1, (2) 



Z 2 L X 1/3 



Air 

where we use the interaction parameter T in the form T = e 2 n 1 ^ i /T. It needs to be mentioned 
that for one-component plasma the damping of oscillations is entirely absent. 

Below, we apply the general concept of stationary level of quasiparticle excitations (here - 
electromagnetic waves) in plasma: damping of electromagnetic waves is compensated for by 
the permanent creation of the new waves in the equilibrium plasma at constant temperature. 

The number of oscillations in volume V with the components of the wave vector k in the 
intervals dk = dk x dk y dk z is equal Vd s k/(2ir) 3 . Then, the number of oscillations with the 
modulus k = \ k | equals Vk 2 dk/2ir 2 . Taking into account the existence of two independent 
polarizations of oscillations and the dispersion equation for the transversal oscillations in 



plasma uo = ^c 2 k 2 + Q 2 , we find the number of oscillations with the frequency u in the 
form 

V ^ijj 2 — Q 2 iijdijj 



7T 2 C 3 



(3) 



Naturally, the Bose distribution by frequencies is valid for the transversal modes in plasmas, 
which are the Bose quasiparticles. Therefore, the spectral distribution for electromagnetic 
waves in plasma has the form 



^^ssSt*""' 1 ' (4) 

where 9{uj — Q p ) denotes the unit-step function which reflects the fact that transversal 
oscillations with u < Q p are absent in plasma. By introducing the dimensionless variable 
q = Tiuo/T and the parameter a = hQ p /T, we can rewrite (jl]) in the form: 



, pf w, VT A q 2 ^q^a~ 2 

dE{q)/dq = ¥^ exp(q)-l 6{q - a) - (5) 
The changes in spectral curves as a function of q for various parameters a are given in the 
Fig. 1. 

In the plasma, the position of maximum in the spectral distribution is a function of the 
parameter a (Fig. 2). 

For a —>■ the maximum value approaches the Planck limit g^ aa; = 2, 822. For the case 
a ^> 1 the position of maximum approaches q x max ~ a + ^ (according to Eq. ([51). Now, we 
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have q m ax > a, for any a, and not only for the asymptotically large values of the parameter 
a. In the classical limit, we find dE{uj) / duj from Eq. (j3J) 

VujT 



dE{uj)/duj = ^ T Ju 2 -Q 2 p 6{u-Q p ). (6) 



Formula (pEJ) is a generalization of the Rayleigh- Jeans classical result (e.g. [6]) for the case 
of fully-ionized plasma. Formally, Eq. (jSJ) can be easily obtained by calculating the energy 
H of electromagnetic field in a transparent medium [7] 

1 f/rfN d[ue{u))] 2 



4, 1<*V^ + < 7 » 
if we use for the frequency-dependent dielectric function e(u) in plasma (the case of a small 
damping, when lme(u) — > 0) the expression e — 1 — fi 2 /u; 2 and take into account the 
expressions for the correlation functions of the electrical (E 2 )^ and of the magnetic (H 2 )^ 
fields 

e(w) (£ 2 ). = (i? )o; = coth f-j . (8) 

It is necessary to stress here that for the ultra-relativistic temperatures, when T > mc 2 
(here m is the mass of the charged particles which determines the dominating contribution 
in the current plasma frequency of the primordial plasma), the spectrum of the transversal 



mode is different from uj = a/^ 2 + c 2 k 2 and depends on the temperature. The spectrum of 

the transversal mode is determined by the dielectric function 

tr -i 7re 2 nc f 2oj / u 2 \ . , ck — u , ] 
e =1 + -777^ --T+ ln l777T-| > ( 9 ) 



cj/cT [ cfc \ c 2 k 2 J ck + uj 
where the charge e and the density n are related to the dominant (for determination of the 
transversal plasma mode) fraction of the plasma, existing at any given moment. In general, 
the expression for the spectrum cannot be written in analytical form. For the region of 
frequencies uj 2 ^> 2ire 2 nc 2 /T, the radiation distribution is close to the Planck distribution, 
because the spectrum approximately is described by the formula uj 2 ~ 2ime 2 c 2 /T + c 2 k 2 ~ 
c 2 k 2 . The asymptotic behavior of the spectrum for a small k has the form 

2 47re 2 nc 2 3 2l2 

uj 2 = h-c 2 r (10) 

T 2 v ' 

Therefore, the gap in the spectrum of radiation distribution and the deviation from the 
Planck distribution are structurally similar to Eq. (J3J). However, in the case of the ultra- 
relativistic temperatures, the deviation appears near the value of uj 2 = ^L2HL < instead 
of uj — flp being valid in the region of the temperatures T <C mc 2 . 
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J diuuj^uj 2 - Vl 2 p ln[l - exp(-hLj/T)]. (11) 



Let's now calculate the thermodynamic properties of the equilibrium electromagnetic 
radiation in plasma (for simplicity we restricted the consideration below by the region of 
the non-relativistic temperatures, where the analysis can be done in a more general form). 
Since the chemical potential is equal to zero, the thermodynamic potential Q coincides with 
the free energy F, which can be written in the form 

8ttVT ro ° 
(2tt) 3 c 3 

We suppose that the quasiparticles in plasma can be in first approximation considered as 
an ideal gas. Then Eq. ( iTTjl can be rewritten as 

F = / dx[a? - MlJT) 2 ] 3 ' 2 -\ . (12) 

3n^h 3 JhQ p /T 1 K p/ n exp(x) - 1 K } 

Therefore, the free energy can be represented as 

F = F <j)(nn p /T), (13) 

where 

1 r°° 1 

4>m P /T) = - / dx[x 2 - m P /Tff 2 (M) 

7T 4 Jm P /T exp(x) - 1 

The function Fq = —V7r 2 T 4 /45(hc) 3 is the free energy for the free photon gas for the case 
when plasma is absent. For small and large values of the ratio hQ p /T, the function can 
be written as 

<p(hn p /T ->o)~i-JL (nn p /T) 2 , (15) 

and 

45^ 

7T 

respectively. The function <f>(a = hQ p /T) is shown on Fig. 3. 
The entropy of radiation S = —dF/dT is expressed as 



<j)(htt p /T > 1) ~ ^^{hVt p /T) 3 ' 2 exp(-htt p /T), (16) 



T 

where the function \ is determined by 



(17) 



X m p /T) = - / dx[x 2 - m P /Tf] 1/2 n r , (18) 
^ JhQp/T exp^xj — i 
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Since x(0) = 5/2ti 2 , the entropy for the case a <C 1 is equal: 

5a 2 



5 — — — 
T 



1 



8tt 2 



(19) 



It means that for a <C 1 the adiabatic expansion of the radiation field is close to the case 
of the ideal photon gas TV" 10-1 , where 70 = 4/3. The small correction 7' to the value of 70 
(7 = 70 + 7') is a function of the system's parameters 7' ~ —5a 2 /2Air 2 lnV. 

For the opposite case, when a ^> 1, taking into account that xi a 3> 1) = 15exp (—a) /it 4 , 
we find the entropy 

4F 45a 2 exp(-a) _ VT^exp {-Ml p /T) 

Let us apply this theory to matter and radiation in the early Universe at the stage when 
the matter was still in the plasma-like state (primordial plasma), but the equilibrium between 
radiation and plasma had already been reached. We will now determine the time-dependent 
plasma frequency Q p (t) = T, a Qp(t), which is the sum of the different time-dependent plasma 
frequencies of the charged particles existing at the moment t. Time dependence is the 
consequence of the expansion of the Universe in time V = V(t) and, the connected with 
this expansion, transformation of the number and species of particles N a (t) in primordial 
plasma. If one denotes, in the model, the moment of the primordial plasma recombination 
as to, and uses the adiabatic relation for radiation field V^~ X T = const with 7 = 4/3 (which 
follows from the conservation of the Planck radiation entropy Sp), then 

_ 4F _ 47TW* 
Sp ~ ~T ~ 45W (21) 

In reality, the recombination process took time which is small in comparison with the typical 
expansion time. Let us suggest, in accordance with the model of Hot Universe, that in the 
moment to, the Planck distribution for radiation was valid, as the density was too low to 
satisfy the inequality hQ p (to)/T(t ) <^ 1. The accepted estimation for the characteristic 
recombination temperature is T(to) ~ 4000K . Let us consider the behavior of parameter 
a(t) = hQ p (t)/T(t) for earlier time t to- Taking into account Eq. (12"IT) we find 

a(t)=aM (^y° =aM (iuy\ (22) 

Therefore, for the model of the early Universe, the parameter a(t) decreases in time and for 
t <C to can be of the same order or bigger than unity. For simplicity we suggest that the 
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quantity of the lightest (a = a ) charged particles is conserved, i.e. N ao (t) = const. The 
more accurate condition for the decrease of the parameter a(t) in time, which takes into 
account the changing composition of the charged particles in primordial plasma, requires 
the increase in time of the ratio f2 p (t)y 1,/3 (t)/^p(to)V Al / 3 (to)- The condition a(t) > 1 is 
attainable if the period of time between the set-up of the time-dependent equilibrium state 
for radiation field in the primordial plasma and its recombination r , is of the same order 
or bigger than r = t\ — t , where t\ is the root of the equation a(ti) = 1. If the condition 
a(t) 3> 1 is fulfilled, then, we have the isentropic process with the entropy determined by 
Eq. ( l20l) . In this case (if the number of the light charged particles is conserved) we find the 
time evolution of the parameter a(t) for t < t' -C t\ (in the region a(t) ^> 1) in the form 

a(t) = a(t') + In (J^j . (23) 

Therefore, the parameter a(t) decreases in time for the region o < 1, as well as for the 
region where a ^> 1. If the state of the plasma and the radiation field at the moment of 
recombination are characterized by the parameter a(t ) <C 1 and this state of the system 
is the result of the adiabatic expansion of another state at the moment t to with the 
parameter a(t) > 1, we arrive at 



T(to) a 2 (t 



o 



4>{a{t)) + ^P- X (a(t)) 



(24) 



Tit) a 2 (t) I 

The adiabat for the radiation field is shown on the Fig. 4. 

It is easy to see that the process of adiabatic expansion of the primordial radiation with 
the parameter a(t) > 1 would lead to the different future scenarios of the Universe evolution, 
when the equilibrium radiation field is described after Planck, and by using the generalized 
(for the case of arbitrary hfl p /T) distributions. 
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Figure 1: The dimensionless spectral distributions dE(q)/dq for a=0 (Planck distribution), and 
a = 1,2,3 (non-Planck equilibrium radiation, the curves are from the top down) as functions of 
the parameter q = Hlo/T. The parameter a equals a = Ml p /T. 
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Figure 2: The positions of maximums q of the equilibrium radiation density in plasma as function 
of the parameter a. 




Figure 3: The function <f> determining the free energy F = Fo0(a). 
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Figure 4: The ratio a(to) as a function of the values of the parameter a(t) for the fixed value 
of the parameter T(to)/T(t) = 0.0001. The upper curve: the arbitrary values of the parameter 
a; calculation based on non-Planck distribution. The lower curve is calculated with the Planck 
distribution extended to the region a > 1, where this distribution is not valid. The curves are 
plotted in the logarithmic scale. 



